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Abstract. We briefly review theoretical calculations for the pseudoscalar-vector meson hyperfine splitting 
with no open flavor and also report a many body field theoretical effort to assess the impact of chiral 
symmetry in the choice of effective potentials for relativistic quark models. Our calculations predict the 
missing rjb meson to have mass near 9400 MeV . The radial excitation rj c (2S) is in agreement with the 
measurements of the Belle and most recently Babar collaborations. 
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Shortly after the discovery of the J/ip it was under- 
stood that a rich spectroscopy of new mesons awaited clas- 
sification. In this task the constituent quark model was a 
useful tool providing a simple periodic table where spec- 
tra and various radiative decays could be correlated with 
the help of a modest number of parameters. In this pic- 
ture vector mesons are a qq pair, in an S or D wave, with 
spins parallel giving total angular momentum J = 1 . Pseu- 
doscalar mesons correspond to the J = ground state 
with S-wave qq pairs spins antialigned. Ignoring the D- 
wave component, the only difference between both sys- 
tems is the relative spin alignment. Any spectroscopic 
mass splitting can conveniently be incorporated in the 
quark model with a term, A<j\ ■ <r 2 that is reminiscent 
of the electron-nucleus spin-spin coupling, hence the name 
"hyperfine" . This was immediately noted by Appelquist et 
al. pQ who predicted a charmonium splitting, AM (J ftp — 
rjc), of about 65 MeV. They extracted the amplitude A 
by estimating the J/ip electron-positron width, r e - e +, to 
be 4 keV. Using the currently accepted value of 5.3 keV, 
the splitting would be about 84 MeV, or about a factor of 
2 smaller than the accepted experimental value of 3097 — 
2980 ~ 120 MeV. The need for a confining potential |2 
was soon understood and calculations (by Appelquist and 
Politzer, and independently Schnitzer [2]) including a con- 
fining strength yielded a larger splitting (40 — 80 MeV) 
than purely Coulombic potentials (15 — 20 MeV). 

In retrospective we see that many of the early models 
utilized scalar confining potentials, which provided a good 
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spin-orbit coupling and radial excitations, but underesti- 
mated the hyperfine splittings. In the early eighties, and 
with the rj c experimental state now known, this splitting 
became a benchmark for new model calculations ^BSE 
U] which now also predicted the corresponding splitting in 
bottomonium. The variation in these predictions is sum- 
marized in Table ^ Subsequently, further progress was 
achieved through improved, renormalized non-relativistic 
perturbative QCD calculations (NRPQCD) (HlEl which 
described bottomonium as a non-relativistic system. How- 
ever, the calculated radii of most bb states are too large 
indicating that a Coulombic description, where the rela- 
tivistic splittings scale linearly with the quark mass, is not 
reliable and that strong interactions still induce important 
corrections at this scale ^U]. Nevertheless, approximate 
ground state descriptions are feasible and useful for ex- 
tracting c and b quark masses. 

Non-perturbative lattice calculations with large error 
bars have also been performed |HU12j for bottomonium 
which yield about half, or less, the hyperfine splitting ex- 
hibited in charmonium. This again indicates the system is 
not fully Coulombic since the splitting is not proportional 
to the quark mass. 

Extending this analysis accurately to the n-p system 
is not currently feasible for either the perturbative or lat- 
tice approaches. Thus one still relics on constituent models 
where the hyperfine splitting has a 1/M 2 dependence on 
the constituent quark mass 0] . This can describe the large 
n-p splitting but not simultaneously the hadron scatter- 
ing phase shifts O n the other hand, we know that 
the pion's mass is very low because of its Goldstone boson 
nature from spontaneous chiral symmetry breaking. Thus 
it is natural to seek a field-theoretical formulation of the 
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Table 1. Predictions for the splitting between vector and pseu- 
doscalar bb mesons. CQM stands for Constituent Quark Model. 
Units are MeV. 



Date, Authors 


Model 


Splitting 


1983 Godfrey & Isgur 


CQM 


60 


1983 McClary & Byars 


CQM 


101 


1985 Igi & Ono 


CQM Coulombic 


60 


1985 Igi & Ono 


CQM log running 


90 


1989 Song 


CQM 


55 


1994 Eichten & Quigg 


CQM Cornell 


141 


1994 Eichten & Quigg 


CQM various 


87/65/64 


1994(98) Davies et al. 


Lattice 


30-50 


1998 Pineda & Yndurain 


NRQCD 


47(20) 


2000 Lengyel et al. 


CQM 


46 


2003 Ebert et al. 


CQM 


60 



quark model which implements chiral symmetry consis- 
tently. Such an approach would predominantly attribute 
the hyperfine splitting in light mesons to chiral symmetry. 
This permits using a more moderate hyperfine potential 
to then decribe the smaller splittings which are exhibited 
in light meson excited states and heavy mesons, both of 
which are not governed by chiral symmetry. Thus we con- 
sider the Hamiltonian (inspired in Coulomb gauge QCD) 



H eff 

T 



= T + V C + V T 



V c 



/ cZx<P T (x)(- ia. ■ V + m 9 /3)!f r (x) 
~ J dxdyp a (x)y(|x-y|K(y) 
Vt = \J dxdyjf (x)j;(y) x 

#(|x-y|) • 



(1) 
(2) 

(3) 



(4) 



Here p a = ^T a <P and J a = ^a.T a W are the quark 
color density and current, respectively. This Hamiltonian 
has been diagonalized previously ^2] f° r Vr = in the 
Bardeen-Cooper-Schrieffer (BCS) approximation for the 
vacuum. These earlier studies of the gap equation deter- 
mined that the dynamical chiral symmetry breaking from 
only a longitudinal potential is relatively small and yields 
a low condensate (pW)o ~ -(100 MeV) 3 . On the op- 
posite limit, calculations for high quark masses using the 
Tamm-Dancoff (TDA) and Random Phase (RPA) approx- 
imations for both harmonic oscillator |15U16| and linear 
potentials ^3 produce almost degenerate pseudoscalar 
and vector meson ground states. They are thus unable 
to describe the charmonium hyperfine splitting although 
the RPA can reproduce the n-p splitting by sufficiently 
lowering the quark mass according to Thouless theorem. 
More recently a study |18j implementing chiral symme- 
try used Vc — and a contact potential for Vr to obtain 
a link with transverse one-gluon exchange, which is sup- 
pressed in our approach by the large gluon mass gap |19| . 
Because that model does not include radial excitations 
or confinement, we have generalized |2()j the treatment 



by employing both a Coulomb instantaneous interaction 
and a transverse hyperfine potential. For the longitudinal 
Coulomb interaction we utilize a potential derived |21| 
from QCD through a BCS truncation of the gluon sector, 
represented in momentum space by 
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1.41 



for p > m g 



V(p) 



12.25 m 193 

for p < m g 



V 



3.93 



(5) 



This is numerically close to the standard Coulomb + lin- 
ear potential. The transverse potential, due to non-explicit 
Lorentz covariance in Coulomb gauge QCD, can be differ- 
ent. Since this term has not been studied theoretically, we 
proceed phenomenologically and choose the same Coulomb 
tail as in Eq. JSJ. ft is then matched at low momentum 
to a Yukawa representing a massive gluon exchange which 
emerges from intermediate hybrid states in the Fock space 
truncation. Thus we take 



U(p) = C(p) for p > nig 
C h 



(6) 



U(p) 



p 2 + mi 



for p < nig 



The constant Ch matches the potential continuously at 
the nig scale. Thus the only free potential parameter is 
ni g which determines simultaneously the strength of the 
confining term and the logarithmic one-loop running of 
both U and V. We adopt m g = 600 MeV and investi- 
gate alternative transverse potentials in a more detailed 
publication pU] . 

Calculating the resulting gap equation at zero quark 
mass we find a sizeable increase of the BCS quark con- 
densate, to —(178 MeV) 3 , which is now closer to the phe- 
nomenologically accepted values (this quantity is sensitive 
to the high energy behaviour of the potential as previously 
noted by Lagae |22)- In the chiral limit the calculated 
pion mass is effectively zero (numerically a fraction of an 
MeV) and the p mass is about 780 MeV. For the vector 
mesons we include coupled S and D wave channels, since 
the Hamiltonian of Eq. contains a tensor interaction. 

Upon increasing the quark mass, the pion mass growes 
rapidly in the RPA whereas the p mass only slowly in- 
creases yielding the hyperfine splitting plotted in Fig. ^ 
for various meson masses. This figure presents our pre- 
liminary results and reflects the success of this approach 
which incorporates chiral symmetry and is simultaneously 
applicable to a wide range of quark masses. 

For the same model parameters, we also predict the 
mass of the missing rjb state, a most important issue in 
hadronic spectroscopy We concur with NRPQCD 

and lattice studies but predict a slightly larger splitting, 
(see below) of about 60 MeV. Subtracting this from the 
T(9460) mass yields 77b (9400). This decreasing hyperfine 
strength trend with increasing quark mass (see Fig.^l in- 
dicates that the potential is not yet scaleless. Note that in 
both PQCD and our approach (see Eq. EJl a hadron scale 



Felipe J. Llanes-Estrada et al.: Chiral Symmetry and Hyperfine Meson Splittings 



3 



Hyperfine splittings vs. average multiplet mass. 
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Fig. 1. RPA ground state hyperfine splittings, M(D~) — 
M(0 ), versus the average multiplet mass (3M(D~) + 
M(0 l "))/4 (circles). The three diamonds represent the ob- 
served TT-p and rj c -J/tp and the NRPQCD rj^-T splittings. 



appears logarithmically in the coupling constant. Also for 
bottomonium there is a small difference between the RPA 
and TDA hyperfine splittings since the TDA r/i, mass is 
about 30 MeV lower than in the RPA. While insignifi- 
cant when compared to the T(9460 M eV) mass, it should 
be accurately included when evaluating a small hyper- 
fine splitting. Non-chiral preserving models, such as those 
based on Schrodinger's equation, will thus underestimate 
the splitting by at least this 30 MeV. Although this is cur- 
rently comparable to the quoted errors in both NRPQCD 
and lattice calculations (20 — 30 MeV), it may become an 
issue in the future. 

Finally, it is noteworthy that our approach naturally 
extends to radial excitations. For the ip(2S) — rj c (2S) split- 
ting we obtain 56 MeV, in agreement with the BELLE 
|2U result. Note added: the BABAR ^FN collaboration re- 
ports a possible detection of the r] c (2S) corresponding to 
a hyperfine splitting o/55(4) MeV. The agreement is en- 
couraging. 
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87ER40365 (A.S.) and DE-FG02-00ER41135, DE-AC05- 
84ER40150 (E.S.). 
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